A PHOTOMETRIC REDSHIFT OF z ~ 9.4 FOR GRB 090429B 
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ABSTRACT 

Gamma-ray bursts (GRBs) serve as powerful probes of the early Universe, with 
their luminous afterglows revealing the locations and physical properties of star form- 
ing galaxies at the highest redshifts, and potentially locating first generation (Popu- 
lation III) stars. Since GRB afterglows have intrinsically very simple spectra, they 
allow robust redshifts from low signal to noise spectroscopy, or photometry. Here we 
present a photometric redshift of z ~ 9.4 for the Swift detected GRB 090429B based on 
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deep observations with Gemini-North, the Very Large Telescope, and the GRB Optical 
and Near-infrared Detector. Assuming an Small Magellanic Cloud dust law (which has 
been found in a majority of GRB sight-lines), the 90% likelihood range for the redshift 
is 9.06 < z < 9.52, although there is a low-probability tail to somewhat lower redshifts. 
Adopting Milky Way or Large Magellanic Cloud dust laws leads to very similar con- 
clusions, while a Maiolino law does allow somewhat lower redshift solutions, but in all 
cases the most likely redshift is found to be z > 7. The non-detection of the host galaxy 
to deep limits (Y(AB) ~ 28, which would correspond roughly to 0.001L* at z = 1) 
in our late time optical and infrared observations with the Hubble Space Telescope, 
strongly supports the extreme redshift origin of GRB090429B, since we would expect 
to have detected any low- z galaxy, even if it were highly dusty. Finally, the energetics of 
GRB 090429B are comparable to those of other GRBs, and suggest that its progenitor 
is not greatly different to those of lower redshift bursts. 

Subject headings: early Universe - galaxies: high-redshifts - gamma-rays bursts: indi- 
vidual (GRB 090429B) - techniques: photometric 



Introduction 



The burst detections and rapid afterglow identifications of the Swift satellite (jGehrels et al 



2009), combined with intensive ground-based follow-up efforts, have confirme d some gamma-ray 



bursts (GRBs) as am ong the most distant objects known in the universe (jTanvir et al.l [2009J; 



Salvaterra et al.ll2009l ). illuminating the conditions of star formation at the earliest epochs. As 
burst detections push toward progressively higher redshifts, the mere existence of GRBs at these 
times will provide important constraints on models of gravitational collapse, galaxy formation, 
and the early generations of stars. At the same time, high- quality spectroscopy of the burst 



after g lows can be expec ted to reveal element abundances (e.g. IStarling et al.l 120051 : iKawai et al 



2006; iBerger et al.l 12005). host galaxy kinematics, and potentially, the H I fraction of the inter- 
galactic medium (IGM), as the process of cosmic reionization unfo lds (e.g. iBarkana &: Loebll2004l : 
Totani et allhood : iTanvir fc Jakobssonll2007l : iMcQuinn et ailbood ). 



GRBs offer some advantages over other techniques for the selection and study of distant galax- 
ies. Most notably, they have unprecedented luminosity, both of the prompt emission, and afterglow 
(e.g. iRacusin et alj|200a : Bloom et al.ll2009l ). enabling them to provide detailed diagnostics of their 
environments, and pinpointing their host galaxies however faint. However, this utility comes at a 
price - GRB afterglows achieve such brightness only fleetingly, and so the time available to obtain 
redshifts and other information for a burst is often very short (normally < 24 hr). In order to 
realize the ambitions of finding bursts at extreme redshift, and efficiently exploiting high-redshift 
GRBs as probes of this early cosmic epoch, it is necessary to devote increasing effort to the rapid 
identification of GRB near- infrared (NIR) afterglows. In addition to workhorse NIR instrumen- 
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tation at large observatories, a growing number of dedicated facilities and instrume nts have been 



comm issioned, with a primary aim of rapidly locating distant GRBs (e.g. PA IRITEL (IBloom et al 



20061 ); Gamma-Ray Burst Optical and Near-Infrared Detector (GROND) (jGreiner et all 12008 )). 
Follow-up spectroscopy of these cand idates has proved several to be at very high redshift (e.g . 



Kawai et al 



2006 



Salvaterra et al 



2009). 



Greiner et al.ll2009h . culminating in GRB 090423 at z « 8.2 (jTanvir et all 1200 



However, in some cases r apid spectroscopy is no t possible, and we m ust fall back on photometric 
redshift measurements (e.g. iJakobsson et al.l 120061 ; lHaislip et al.l 120061 ) . Here again, GRBs offer 
some advantages over galaxies for the application of such techniques. First, there is little intrinsic 
variation in the spectral shape of an afterglow - it can be modeled simply as a power-law plus host 
galaxy dust extinction Ly-a absorption. This is in contrast to the diverse spectra of galaxies, which 
can have contributions from young/old populations (or a mixture), dust in complex configurations 
, exhibit intrinsic curvature, Balmer breaks etc, none of which are a concern for GRB afterglows. 
Second, the identity of a GRB afterglow is unambiguous from its fading, and thus there is no chance 
of mistaking a GRB afterglow with e.g. a Galactic L or T dwarf, which can also confuse high-z 
galaxy s earches. It has been shown that GRB photometric redshifts are generally robust for these 
reasons (jKriihler et al.ll201ll ). Indeed, while the fundamental accuracy is limited by the bandwidths 



and bands used, GRBs are much less subject to the "catastrophic" failure of photometric redshift 
determination, that can impact individual galaxy measurements. 

In this paper we discuss the discovery and multi- wavelength follow-up of GRB090429B. The 
afterglow was not visible in deep early optical imaging, but was found in deep IR observations 
starting ~ 2.5 hr after the burst. While spectroscopic observations were curtailed by poor weather 
conditions, our photometry does allow us to construct a spectral energy distribution (SED) for the 
burst, and to infer a photometric redshift of z ~ 9.4 , making GRB 090429B one of the most distant 
objects known to date. 

The paper is structured as follows: in Section[2]we present our full dataset on GRB 090429B and 
the uncertainties of our photometric measurements; in Section [3] we derive our photometric redshift, 
supplemented with deep host observations. Finally, in Section H] we summarize our conclusions, 
highlighting the importance rapid-response NIR imaging and spectroscopic capability on large 
telescopes for the study of the early universe using GRBs. Throughout this paper we assume 



ACDM cosmology with H = 72 km s -1 Mpc 



-l 



n 



M 



0.27, f?A = 0.73, and use a standard 



nomenclature to describe the variation of the afterglow flux density as F oc t a v ° . 
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2. Observations and Analysis 



2.1. Swift Observations 



The Burst Alert Telescope (BAT; iBarthelmv et al.ll2005l ) aboard the Swift satellite triggered 
on GRB 090429B at T = 05:30:03 UT. The 15-350 keV light curve is composed of three distinct 
peaks with a total duration T90 = 5.5 s, and the time-integrated spectrum can be fitted by a 
single power law with an exponential cut off. The derived total fluence in the 15-150 keV band 
is 3.1 x 10 _7 erg cm~ 2 , with E pea ^ = 49keV. This peak energy is among the few detected by Swift 
within the BAT bandpass. 

After 106 s, the narrow-field inst ruments began their standard burst-response observation se- 
quence. The X-ray Telescope (XRT: iBurrows et al.ll2005l ) identified an uncataloged fading source 
at RA(J2000)=14 h 02 m 40!10, D ec(J2000)=+32°10 / 1 4 / . / 6; no optical/UV counterpart was seen in the 



UV-Optical Telescope (UVOT: iRoming et al.ll2005l ) data. The X-ray data has been characterized 
using standa rd routines in Hea soft, Xspec, and QDP, with the light curve fitting process as 
described in lEvans et al.l (|2009l ). For some analys es, we have used t he au tomatic data products 
produced by the UK Swift Science Data Centre (|Evans et al.1 120071 . l2009l ) . Our presentation of 
parameters derived from the Swift data follows the convention of quoting errors at the 90% level. 
The time-averaged 0.3-10 keV X-ray spectrum from 97 to 29893 s after the burst is best fit by a 
power-law with photon spectral index Tx = 2.01^Q24 and with a total absorption column density 



10 20 

cm 2 , mildly (2.7<r) in excess of the Galactic absorption of 1.2 x 10 20 cm 



of iVn 

we discuss the possible significance of this finding in Section 

The X-ray light curve, given in Table CD and illustrated in Figure CO is adequately fit by a 
combination of brightening and fading temporal power laws: initially, the X-ray flux rises with 
temporal index axi = — 0.96lto 52' referenced to the burst time; following the peak time Tx = 



589 



+146 
-80 



s, the light curve then breaks to a power-law decay with 0x2 = 1-20 



+0.08 
-0.07- 



2.2. Optical and Near-IR Observations 

Basic reduction steps for all optical and NIR photometry were performed using IRAF soft- 
war^E Photometric analysis used both IRAF and the Starlink GAIA software, as well as our own 
custom scripts. Errors in the sky subtraction step are estimated from multiple apertures of size 
equal to that of the source aperture, placed around the field of the GRB. 

Optical images were calibr ated using field stars from the Sloan Digital Sky Survey (SDSS) 
Data-Release 7 (DR7) catalog (jAbazajian et al.ll2009l ) in the and NIR images were provisionally 



1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation 
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calibrated directly to the Two Micron All Sky Survey (2MASS) catalog, but subsequently refined 
as described below. Detections and limits on the brightness of any associated source are presented 
in Table H 



2.2.1. ES02.2m/GROND Observations 



The GROND (jGreiner et al.l 120081 ) observed the field of GRB090429B simultaneously in its 
(dichroic and filter d efined) g'r'i'z' JHK S filter set beginning 14 minutes after the Swift discovery 
(jOlivares et al.ll2009i ). No source was detected at the X-ray afterglow position in any of the seven 
bands: the limits being shallower than usual due to the high airmass for this (northern) field. 
Nonetheless, the implied X-ray to optical spectral slope of /?ox < 0.1 implied suppression of the 
optical flux relat i ve to the X-ray, rendering GRB 090 429B a "dark" burst under the definitions of 



Jakobsson et all (|2004l ) and Ivan der Horst et al.1 (|2009i ). 



2.2.2. VLT Observations 



Deep R and z-band observations were made with the VLT/FORS-2 camera at ~60 minutes 
post-burst. Once again no optical source was visible at the position of the X-ray afterglow, con - 
firming that it was unusually dark, and thus a good candidate high-z GRB (jD'Avanzo et al.ll2009l ). 



2.2.3. Gemini-North Observations 



Beginning roughly 2.5 hr after the burst trigger, we carried out a series of observations 
from Gemini-N orth. We gather ed optical i' z' imaging with the Gemini Multi-Object Spectro- 
graph (GMOS; Irlook et al.l 120040 and NIR JHK imaging with the Near-Infared Imager (NIRI; 
Hodapp et al.ll2003l ). GMOS observations consisted of five exposures of 3 minutes each, per fil- 
ter; NIRI observations consisted of eight dithered positions of 60s each. The Gemini GMOS and 
NIRI packages under the IRAF environment were used to sky-subtract, align, and combine the im- 
ages. The NIRI images were also corrected for the small non-linearity effect seen in the detector:^. 
Photometry was performed relative to SDSS stars for the GMOS data, and relative to secondary 
calibrators from GROND for the NIRI data (see Section I2.3() . Our photometry is presented in 
Table [2j 

While no optical counterpart was present in our i! or z' images, we did identify a source within 
the X-ray localization in our NIR observations. The position of the source was RA=14 h 02 m 40.10 s , 
Dec =+32°10 / 14'.'20. Following this discovery we attempted spectroscopic observations from Gemini- 



2 http://www.gemini.edu/sciops/instruments/niri/data-format-and-reduction/detector-linearization 
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North, however, increasing summit winds forced the closure of the telescope and meant that these 
were aborted with <10 minutes of useable exposure time and no trace is visible in the observations. 
We obtained a second epoch of K-band observations on April 30 UT revealed a clear fading of 
~1.2mag of the identified source, confirming its transient nature (and corresponding to a power- 
law index of ax = 0.53 ± 0.10, shallower than the X-ray decay at that time). Figure [2] presents 
our Gemini imaging data, while the lower panel of Figure [T] shows our optical/NIR lightcurve. The 
resulting SED, from X-ray to IR is shown in Figure El 

No evidence of a host galaxy is present in our images. A deep r'-band image of the field, 
taken again with GMOS under good conditions (0.4" seeing) at 14 days after the GRB, is shown in 
Figure HI This allows us to place a 3a upper limit on the host galaxy apparent magnitude of r' > 
27.07 mag. We also note in these images the presence of a massive elliptical galaxy, offset roughly 
45" fro m the GRB locati on. This galaxy has absolute magnitude M r ~ —21.6 and Mk ~ —24.5 
(~ L*; iJones et al. IBS)- It appears to be the central galaxy of a modest cluster at z = 0.07fl 



It is likely that this foreground structure provides some lensing boost to the observed flux of the 
burst, although the relatively large impact parameter suggests it will not be a major factor. 



2.2.4- HST Observations 

We obtained late time observations of the field of GRB 090429B with the Hubble Space Tele- 
scope (HST). These were taken after the afterglow had faded, and had the goal of finding or 
constraining the host galaxy. We used both the Advanced Camera for Surveys (ACS) and the 
Wide Field Camera 3 (WFC3). Observations were obtained in F606W (broad V - R), F105W 
(broad Y — Z), and F160W (H): a log is given in Table El The data were reduced in the standard 
fashion using multidrizzle and the HST archive "on-the-fly" calibration. All the images were 
drizzled to a common pixel scale of 0"05 pixel -1 . We ascertained the location of the burst on 
the HST images via relative astrometry between our first epoch i^-band observations, and those 
obtained with HST. Doing so we used a total of 11 and 10 sources in common to each frame, for 
ACS and WFC3 respectively. The resulting astrometric accuracy is C/08 (F606W), / /07 (F105W) 
and C/06 (F160W) respectively. At the location of the afterglow we see no obvious host galaxy 
candidates in any of the images. To quantify the depths of these images we estimate the sky 
variance from a large number of background apertures (~50) placed in the field around the target 
position, avoiding visible sources. We then measure the resultant flux at the target position in an 
aperture of Cf /4 diameter, consiste n t with the approaches of many groups in searching for high-z 



galaxies (e.g. iBouwens et al.l l2010i . l201ll ). Our fluxes are shown in Table El In addition to the 
measured fluxes we also show the effective AB-magnitude limits at these locations, which are equal 
to the measured flux density + 3cr, with an additional aperture correction to account for light 
missing within our small measurement apertures. These corrections are small for ACS (0.18 mag 



3 Redshift and r magnitude of galaxy obtained from the SDSS DR7 database; K magnitude from 2MASS. 



for F606W, ISirianni et al.l 120051 ) . but larger for the WFC3 images (0.31 and 0.54 magnitudes for 
F 105 Wand F160W respectively 



2.2.5. UKIRT Observations 

We obtained observations with the United Kingdom Infra Red Telescope (UKIRT), Wide- 
Field Camera (WFCAM), beginning April 29 at 09:18 UT, roughly 4 hr post-burst. Only a limited 
number of exposures were possible due to high wind keeping the telescope shut much of the night. 
These observations were not deep enough to reveal the afterglow, however because the large field of 
view (13.6 arcmin square for each chip) includes many bright 2MASS stars these images allowed us 
to precisely determine the magnitudes of fainter stars, which was crucial for calibrating the NIRI 
images (see Section [2.3p . Pipeline reductions were performed by the Cambridge Astronomical 
Survey Unit (CASU@)~^ 



2.3. Precise Photometric Calibration and Uncertainties 

Since our photometric redshift analysis will depend critically on the accuracy of our photom- 
etry, we took particular care in both the calibration and estimates of photometric uncertainties. 

The Gemini-North/NIRI detections are crucial, but also difficult to analyze since the field of 
view is small (2 arcmin on a side) and there is only one 2MASS star (namely star B in Table [U 
and Figure 2]) that is in all the sub-exposures of the nine-point dither pattern. There is another 
2MASS star (A in Table H|) which appears on two of the sub-exposures, and we used this as a 
double check on the derived photometry. Both these stars are toward the faint end of the 2MASS 
catalog and have relatively large photometric uncertainties. To overcome this we used the wide-field 
UKIRT /WFCAM and ES02.2m/GROND JHK images (both of which were obtained close in time 
to the Gemini observations), which were very precisely calibrated using many bright 2MASS stars, 
to obtain more accurate magnitudes for these reference stars. The two independent determinations 
were consistent with each other within their respective calculated errors (typically 0.01-0.02 mag), 
and we therefore formed a weighted average to obtain our best estimates of the NIR magnitudes, 
as shown in Table [2j 

Magnitudes for the afterglow were measured relative to star B, although this procedure was 
further complicated by the fact that the point-spread-function (PSF) was found to change across 
the frame resulting in the core of the reference star becoming noticeably extended when it was close 
to the southern edge of the detector, as it was in some sub-exposures. This precluded small aperture 



4 http: / /www.stsci.edu/hst / wfc3/phot_zp_lbn 
5 http://casu. ast.cam.ac.uk/ 
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(5 pixel radius, « C//6 ) photometry for these exposures, so in such cases we used a fainter star (C 
in Table [2]) closer to the GRB position as a secondary reference, having determined its magnitude 
relative to star B using those frames where it was not near the edge. We note that a small aperture 
was required to maximize the signal to noise for the afterglow, and that profile-fitting photometry 
was deemed inappropriate due to the small number of bright stars available to define the PSF. 

The magnitudes (and errors) for the afterglow in each band were then determined from an 
error weighted mean of the differen t sub-exposu res. Finally, we converted to flux density using a 
recent NIR spectrum of Vega (see iBohlinl 120071 ) which resulted i n values that are 2% -3% higher 
for our passbands than found using the conversion in Table 7 of iHewett et al.l (|2006l ). These are 
the flux densities reported in Table [2j although we note that when we come to the photometric 
redshift analysis (below) we fit in counts rather than flux, to allow for the different spectral shapes 
of the afterglow and comparison star. 

Since the optical observations provided only upper limits the overall fit is not strongly sensitive 
to the precision of the optical photometry. However, in this case, the field of GRB 090429B lies 
fortuitously within the SDSS survey area, and our most constraining optical limits (from GMOS) 
are obtained in the same filter set. This allows a precise photometric calibration of these images. 
For our Very Large T elescope (VLT) obs ervations we calibrate the field using SDSS observations 
and the transforms of I Jester et al.l ([2005). These latter values were confirmed as reasonable using 
archival zeropoints. 



3. Results and Discussion 

3.1. Temporal Behavior 

Since our observations were taken over several hours, temporal variations in the afterglow 
luminosity could effect our analysis. Unfortunately we have rather little handle on the variability 
at optical/IR wavelengths. Although most GRBs begin power-law decline in luminosit y fairly early. 



in so me cases flat or even increasing luminosity can be seen for a period of time (e.g. iRvkoff et al 



2009). A rapidly fading afterglow (similar to those commonly observed) would imply even more 
stringent upper limits in our (earlier) blue band filters since the extrapolation to a common time 
would yield a more extreme limit on color index in e.g. z — H. Alternatively, the more unusual 
case of a rapidly rising afterglow would yield somewhat weaker constraints since the non-detections 
in the bluer bands could be ascribed to the brightening of the afterglow in the time frame between 
the optical and IR observations. However, this is countered by the fact that such a rising afterglow 
would also imply an even bluer H — K color, more difficult to attain with extinction. In this 
scenario, the rising of the afterglow may offer some support for a high-redshift scenario, since the 
time dilation at z ~ 9 would result in a forward shock which takes a factor of ~ 10 longer to reach 
maximum than at z ~ 0. 
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In an attempt to constrain the temporal slope of the optical afterglow we first perform pho- 
tometry on the individual NIRI frames. We find no statistically significant variation over the ~ 10 
minute time frame of these observations, implying that the afterglow is not varying especially 
rapidly. Second, we utilize acquisition images taken prior to the aborted NIRI spectroscopic obser- 
vations. These suggest a minor brightening of the afterglow between 13000 s and 17000 s after the 
burst (0.3 ±0.2 mag, corresponding to a ~ —1.0 ±0.7). However, these observations were obtained 
at a single dither position, and contained substantial persistence. Hence, we can accurately remove 
neither sky nor dark current and the resulting observations contain large variations in the sky on 
relatively short length scales. Thus, we caution against their use for detailed photometric work, 
aside from noting that suggest that the afterglow is neither rising, nor falling at an unusually rapid 
rate. We gain a much better handle on the decay between the first and second night observations, 
which gives a ~ 0.6 ± 0.1, but of course this may not apply during the first few hours post-burst. 

On balance, then, we favor photometric fitting in which the observed magnitudes are assumed 
to be constant over the period of the early observations (i.e. we assume a = 0). This is consistent 
with the relatively flat X-ray behavior between 1 and 3 hr post-burst (Figure H]). For completeness, 
we have included a single power-law temporal decay as a possible parameter within these fits, 
and confirm that for any reasonable slope — 1 < a < 1 our results are broadly insensitive to the 
assumed value of a (see below). To avoid extrapolating over too wide a range of times, and to 
counter against unusual afterglow behavior (which is normally most notable in the first few hundred 
seconds in the rest frame) we also fit only data taken after 4000 s. The inclusion of earlier data 
would further strengthen our results if we assumed the afterglow to be decaying, but would make 
minimal difference to the fit if we assumed a flat or rising light curve (since the early limits are 
shallower than those obtained at later times). 



3.2. Photometric Redshift Analysis 

Here we attempt to derive the redshift of GRB 090429B via our broadband photometry of its 
afterglow. The absence of any detections within the optical window, if interpreted as the signature 
of high redshift, immediately implies z > 6.3. Similarly, if we interpret the red J — H color 
of the afterglow as indicative of a Lyman-break lying within the J-band, the inferred redshift is 
8.0 < z < 10.5. 

To obtain stronger constraints on the redshift of GRB 090429B we performed the following 
analysis. We considered just the seven deepest observations, namely those obtained at the VLT 
and Gemini-North on the first night with filters redder than 6000 A. We assume initially there is no 
temporal variation over the course of observations, although including plausible variability within 
our fits also confirms that our results are broadly insensitive t o this assumption (se e below). After 



correcting for Galactic foreground extinction (Eb~v = 0.015; ISchlegel et al.lll998l ) we fitted these 
flux density measurements with a grid of simple models for the SED of the afterglow. The errors 
are likely to be Gaussian distributed, to good approximation, since the uncertainties are dominated 
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by background subtraction (although we also included zero-point calibration uncertainties in the 
modelling), and therefore used minimum-x 2 fitting. 

Specifically, the model was a simple power law, with the spectral index, /3q, and overall nor- 
malization as free parameters. The grid of models spanned a range in redshift of < z < 12 and 
rest-frame F-band extinction of < Ay < 12. Beyond z ~ 7 we are effectively fitting only three 
data points so there exists a degeneracy between extinction (Ay) and /?o- We therefore include a 
weak prior for the probability distribution of the value of f3o, (see Figure [5]) which is modeled as a 
log normal with a peak likelihood at f3o = 0.5 and a width such that the relative likelihood is 50% 
of the maximum from about 0.3 < /?o < 0.85. This is physically motiv ated since it allows y alues of 
flo over a broad range, comparable to the range usually observed (e.g.. ISchulze et al.ll201ll ). but in 
particular prefers /3o = /?x — 0.5 sa 5, as would be e xpected if there was a cooling break between 
the X-ray and optical regimes (e.g., ISari et al.lll998l ). The plausibility of such a cooling break is 
clear from Figure [3j while relaxing this assumption does not lead to any viable fits at low redshift. 
Added to this was absorption due to neutral hydrogen in the IGM (jMadaulll995l ) , (neutral hydrogen 



in the host was taken to have a typical column density of 10 cm" , although the results are in- 
sensitive to the exact number assumed), and extinction due to dust. We experimented with several 
dust la ws, from th e Milky Way (MW), Large Mage llanic Cloud (LMC) an d Small Magellanic Cloud 
(SMC) (|Peilll992h . as well as the extinction law of lMaiolino et all (j200J). 



We also impose a weak prior o n the intrinsic luminosity of the optical afterglow (Figure [5]). 
Studies such as that of iKann et al.l (|2010l ). indicate that there is an upper envelope to the (broad) 
distribution of GRB optical afterglow luminosities. We therefore apply a prior which is flat below 
this envelope and cuts off exponentially at brighter luminosities, although the cut-off is slow enough 
to allow a reasonable probability that the luminosity could be somewhat higher. 

In fact, rather than fitting directly to the flux densities, we integrated our model spectra over 
response functions and compared the counts obtained by integrating an approximate spectrum of 
the comparison star. The response functions were obtained from the measured filter transmission 
curves, multiplied by a typical atmospheric absorption curve generated by ATRA^. Going to 
these lengths effectively corrects for the small difference in the SED shape of the afterglow from 
the reference star, although again the conclusions are not greatly affected. In Figure [6] we show 
the photometric data points and the best-fit model for the afterglow spectrum assuming that the 
afterglow did not evolve temporally during the first 3 hr (see below) and that the dust is similar 
to that in the SMC, which has frequentl y been found to be a goo d approximation to the dust laws 
along many other GRB sight lines (e.g.. ISchadv et alj|2007l . l2010l ). We also show the best fit low-z 
model (as it happens z ~ 0), which is formally ruled out at high confidence. In Figure [7] we plot 
contours of \ 2 ° ver a grid of models spanning a range of redshift and rest-frame l^-band extinction, 
Ay. The red-cross shows the best fitting model, which has z = 9.36 and extinction Ay = 0.1, 
although the 99% confidence contour runs as low as z ~ 7.7 if there is a modest amount of dust 



3 http://atran. sofia.usra.edu/cgi-bin/atran/atran.cgi 
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(vest-Ay ~ 0.5) in the host. Marginalizing the likelihood (which we define C oc exp(— x 2 /2)) over 
Ay (assuming a flat prior) indicates a 90% likelihood range of 9.02 < z < 9.50. There is no solution 
at lower redshifts (z < 7) which is not ruled out at ^>99.9% level: and the best fit at low redshift 
(z ~ as it happens, as shown by the blue cross) requires very high extinction of Ay ss 10. 

In Figure [8] we show similar likelihood contours for fits spanning a broader range of models 
with different prior assumptions for the temporal power-law decline index a, and commonly used 
dust laws. Changing a to ±1 makes rather little difference, and in any case, as discussed above, 
there is evidence to suggest the luminosity was not changing even as rapidly as thi s. Varyin g the 
dust law does have m ore effect, largely due to the 2175 A feature in the MW, LMC (|Peilll992l ). and 
Maiolino et al.l (|2004l ) laws producing the blue H — K color ev e n at s lightly lower redshifts, although 
in most cases the best fit remains z > 9. The IMaiolino et al.1 (|2004l ) dust law was determined from 
observations of a quasar at z = 6.2 and is argued to consistent with dust produced largely from 
early supernovae (note that this law is only defined up to ~ 3200 A in the rest frame, and we 
therefore graft it to the SMC law at this point). This case is interesting as it does allow redshifts 
as low as z ~ 6.5 at 99% confidence, althoug h to date, only G RB 071025, with a photo- z ~ 5, has 
shown evidence of requiring such a dust law ( Per lev et al.ll2010l ). 



Finally, in Figure[9]we show the likelihood as a function of redshift for the SMC dust-law models 
having marginalized over both a (assumed a flat prior between -1 and +1) and Ay (assumed a 
flat prior between and 12). The maximum likelihood is at z = 9.38, and 90% of the likelihood is 
between 9.06 < z < 9.52. 



3.3. Implications of the Absence of a Host Galaxy 



Our late time date taken with Gemini and HST are potentially extremely valuable, since we 
can use the absence of any host galaxy candidates to assess the plausibility of any lower- z solutions 
to our photometric redshifts (the HST images are shown in Figure [10]) . The detection of a host 
galaxy in the optical was used, for example, to show that GRB060923A was z < 3 despite its 
afterglow being a ET-band drop-out (jTanvir et al.ll2008l ). In the case of GRB 090429B, the possible 
low-redshift scenarios seem to be those with z < 1 and high dust extinction Ay ~ 10 (although 
we emphasize that such models remain formally ruled out). The limits these data provide on this 
are shown graphically in Figure [TT| where we plot the absolute inferred magnitude of the host 
galaxy in the observed V, Y and H bands as a function of redshift. For completeness we cut each 
line at the point where 1216Ax(l + z) passes the central wavelength of the band. At z = 0.1 
close to the minimum of our lower redshift solution we obtain inferred absolute magnitude limits 
of My > —10.6, My > —9.9, Mh > —10.5, these exceptionally deep limits are comparable to the 
luminosities of bright globular clusters, and signifi cantly fainter than a ny known GRB or supernova 
host galaxy, indeed they place limits of < 10 -4 L* ( Blanton et al. 2003). E ven at z = 1 the observed 
Y-band limits would imply M B > -15.1, or < 0.001L* (|Rvan et al.ll2007h . 
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In this regard it is worth noting that the lower redshift solutions are only viable in cases where 
the host galaxy extinction is high, whereas such faint galaxies typically have low metallicity, and 
little dust, and it is therefore extremely unlikely that one could create the extinction (Ay ~ 10) 
necessary to explain GRB 090429B. Any z < 3 solution would require the host of GRB 090429B to 
be fainter than the large majority of GRB hosts currently known. Furthermore, our wide wavelength 
coverage would also allow us to uncover any very red dusty host galaxies, wh ich would provide the 



necessary extinct ion, but would be missing from optical only searches (e.g., ISvensson et al. 



2010 



Levan et alJl200d ) 



This offers strong support for our high-z model, In these cases only the Fl 60 W observation 
yields potential information as to the magnitude of the host galaxy. The inferred 1500 A absolute 
AB magnitude at z ~ 9.4 is —19.95. This lies roughly in the middle of the observed absolute 
magnitude distributio n of z ~ 8 candidate galaxies found in deep HST ACS and WFC3 imaging 
( Bouwens et al. H2010I ). and thus the non-detection of a galaxy at z ~ 9.4 is not unexpected in 



observations of this depth. 



3.4. Other Indicators of High Redshift 



In addition to the above discussion there are additional lines of evidence which offer support 
for the high-z interpretation of GRB 090429B. The best z < 5 solution is actually at z ~ 
(although it remains a very poor fit to the available optical data) and requires Ay ~ 10, for an 
SMC extinction law, which would correspond to a foreground Nr ~ 10 23 cm -2 , nearly two orders 
of magnitude larger than is observed. While the dust-to-gas ratios observed through the MW can 
show moderately large variations, such a large offset would be unheard of, particularly in GRB 
afterglows where typicall y the ratio of dust e xtinction to X-ray determined gas column is actually 
less than is seen locally (jSchady et al.ll2010l ). Hence the observed X-ray spectrum seems to rule 
out any low redshift {z < 1), high extinction scenario. This is illustrated in Figure \7\ which shows 
that the best low-redshift solutions are well above the contours of Ay inferred from the excess Nn 
assuming typical GRB dust-to- gas ratios. 

A second line of evidence comes from observed high-energy correlations seen in many GRBs. 
In parti cular the relation b etween the peak energy of the vF v spectrum (E p ) and the burst isotropic 
energy (jAmati et al.ll2008l ). Although this relation has a significant scatter, it can also be used to 
place some constraints on the burst redshift under the assumption that all long bursts should follow 
the relation. For GRB 090429B, the burst is only consistent with this relation at better than 3<r 
if z > 1, implying that it also disfavors very low redshift models for the origin of GRB 090429B. 
This result couples with the lag-luminosity relation, which would impl y an isotropic lumino sity of 



Lj so ~ 10 53 erg s , similar to other long GRBs in the "silver sample" (jUkwatta et al.ll2010l ). 
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3.5. Rest-frame Properties 



The observed fluence of GRB 090429B is 3.1 x 10 erg cm , comparable to that observed for 
GRB 090423, and implies an isotropic energy release in the 15 — 150 keV band of E- lso = 3.5 x 10 52 erg 
at z ~ 9.4 . Its absolute X-ray (rest-frame ~3-100keV) brightness at ~ 1000s of ~ 2 x 10 49 ergs -1 , 
and i^-band luminosity of M\/(i+ z ) of ~ — 26.2(AB) are also similar both to GRB 090423 and the 
bulk of the long GRB population. 

The rest-frame duration at first sight seems surprisingly short, with Tqq/(1 + z) ~ 0.5 s. Inter- 
estingly, three of the four highest redsh ift GRBs discovered prio r to this one have also had rather 



low values of Tqq/(1 + z) a round 1-4 s (jRuiz-Velasco et al 



2009 



2007 



Greiner et al. 



2009 



Tanvir et al 



Salvaterra et al.ll2009i ). A possible explanation for this tendency may be that for high redshift 



sources the BAT is observing at rest frame MeV energies, where the light curves tend to be more 
rapidly variable and shorter duration than at lower energies, and thus it is plausible that just a 
single peak of emission is detected rising up above the noise in these cases. 



An analysis obtained using the technique discussed in lUkwatta et al.l (|2010l ) reveals that this 
GRB presents a small value of lag31/(l + z) = 58 ± 27ms between the low-ener gy (15-25 keV) 



and th e high-energy (50-100 keV) channels. The rest-frame correlations found by lUkwatta et al 
(|2010l ) indicate this corresponds to a peak isotropic luminosity L lso > 10 52 ergs _1 , consistent with 
the observed fluence and short duration. 

Another interesting issue is that of the absorption inferred from the X-ray spectrum. Although 
the measurement is not highly significant, if taken at face value, the rest-frame column density is 



10 23 

cm 2 (90% confidence range). This would be already high compared to 



most other Swift observed GRBs, and would be higher still if, as is very likely, the metallicity is 
substantially less than Solar (which by convention is often assumed in calculating Nu). Such a 
high-column density is no doubt surprisi ng, although a similar value was found for GRB 090423 
(jTanvir et alJ 12009 ; ISalvaterra et al.ll2009i ). As in that case, it also raises the question of whether 
a high gas column would be compatible with the low extinction indicated by the NIR afterglow. 
GRBs are expected to b e able to destroy dust to fairly large distances from their birth sites (e.g., 
Waxman k Drainel I2OO0I : iFruchter etallboOll ). but the good fit for many afterglow SEDs with a 
SMC extinction law suggests they have not generally bee n highly modified b y dust destruction, 
which would tend to produce "gray" extinction laws (e.g. ISchadv et al.ll2010l ). In any event, the 
large error bar on this measurement and potential systematic uncertainties in calibrating the soft 
response of the XRT, taken together with the wide range of dust-to-gas ratios seen to other GRB 
sight-lines, makes the significance of this conflict hard to assess at the present time. 

On balance, it seems that in most respects the general properties of GRB 090429B do not stand 
far apart from the population of long-duration GRBs, even at the inferred redshift of z ~ 9.4 . 
In particular it shows no evidence that its progenitor is distinct from those of GRBs seen in the 
more local universe. This is of particular importance at z ~ 9.4 , since it is close to the redshift 
where WMAP observations imply the bulk of reionization of the universe occurred (z = 10.6 ± 1.2, 



Komatsu et al.ll201ll ) . This reionization process is likely to have been driven by the first generations 



of star formation, including Population III whose pristine H+He composition is expected to lead 
to generally more massive stars. It has been proposed that a consequence of t his could be that 
Popu lation III stars produce particularly long duration, and energetic GRBs (e.g. iMeszaros fe Rees 



201Q ). This is clearly not the case for GRB 090429B, and hence we conclude that its progenitor 



was more likely to be a high-mass second generation (Population II) star. 



Conclusions 



We have presented our discovery and multi- wavelength observations of GRB 090429B and its 
NIR afterglow, and a deep late-time search for its host galaxy. The afterglow exhibited a strong 
spectral break in the J band, which coupled with the non-detection in the optical, and relatively 
blue H — K color, allows us to derive a best-fit photometric redshift of z ~ 9.4 . It is, of course, 
important to look carefully at the evidence against a lower redshift origin, since we know that Swift 
GRBs exist in much greater number at z < 4 than above it. Our afterglow photometry allows us to 
exclude all z < 6 solutions with high confidence. A low redshift (z < 1) is also effectively ruled out 
by our HST observations which would easily locate such dusty galaxies in either the optical or IR 
at z < 1, and also the relatively modest excess iVn which would not be consistent with a very high 
dust column. This immediately implies that GRB 090429B is one of the most distant objects yet 
discovered. The maximum-likelihood solution, with our preferred assumption of an SMC extinction 
law, is z = 9.38 with a 90% likelihood range of 9.06 < z < 9.52. The conclusion s do n ot depend 
sensitively on the priors adopted for other parameters, although a lMaiolino et al.l (|2004l ) dust law 
would favor somewhat lower redshift, but still z > 7, at the expense of requiring fairly significant 
extinction (up t o rest frame Ay ~ 2). Since all z > 6 bursts observed to date are consistent with 
having Ay = (jZafar et alJl2010l . l201ll ). this suggests that z ~ 9.4 provides a good estimate of the 



redshift of GRB 090429B. 

Our campaign shows again how rapid-response multiband NIR observations play a crucial role 
in identifying candidate extreme-redshift afterglows. However, it also highlights the need for even 
more rapid observations and decisions to maximize the likelihood that spectroscopic observations 
can be successfully obtained. In the future, additional dedica ted ground-based optical/NIR multi- 
band imagers such as GROND and RATIR (jFarah et al.ll2010l ) can be expe cted to feed further suc h 
candidates directly t o NIR spectrograph s including X-Shooter on th e VLT (ID'Odorico et al.ll2004h . 
FIRE on Magellan (jsimcoe et al.l hood ), and GNIRS on Gemini (jElias et al.l hood ); ultimately, 
such prompt spectroscopy of extreme-redshift candidates will not only resolve the nature of these 
events, but quite likely succeed in realizing the extraordinary promise of GRBs as probes of the 
extreme-redshift universe. 



The Gemini data, acquired under the program ID GN-2009A-Q-26, are based on observations 
obtained at the Gemini Observatory, which is operated by the Association of Universities for Re- 
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Table 1. X-ray Observations 



T — Tn (s) 


"Plnv P)pnsit,v ( //.Tv ) 

-L 1UA J-yv^llOlU y \ fUitJ y I 


Prvor ( //.Tv ) 


158 


0.53 


0.12 


291 


0.80 


0.18 


368 


1.18 


0.27 


452 


1.19 


0.27 


561 


1.69 


0.36 


651 


1.37 


0.31 


708 


1.35 


0.30 


768 


1.18 


0.27 


838 


1.00 


0.22 


927 


1.31 


0.21 


4618 


0.099 


0.026 


5636 


0.111 


0.029 


6417 


0.126 


0.033 


10316 


0.089 


0.023 


10982 


0.087 


0.023 


15921 


0.0269 


0.0053 


120249 


0.0021 


0.0006 



Note. - X-ray observations obtained by the 
XRT instrument onboard the Swift satellite. The 
flux density is calculated at 2 keV, while the con- 
version factor from flux to Jy is 8.87 x 10 erg 
cm~ 2 s" 1 Jy -1 (0.3-10 keV). The counts to flux 
conversion factor is 2.984 x 10~ n erg cm -2 s _1 
count -1 . Uncertainties are la. For the best fit 
X-ray spectrum parameters see Section [2j 
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Table 2. Log of Ground-based Optical/NIR Observations 



T - T (s) 


Magnitude 


Flux Density (//Jy) 


Filter 


Telescope 


990 


> 23.08 




g' 


GROND 


990 


> 22.86 




f 

r 


GROND 


990 


> 22.03 




a 


GROND 


990 


> 21.87 




z' 


GROND 


990 


> 21.06 




j 


GROND 


990 


> 20.50 




H 


GROND 


990 


> 19.84 




K 


GROND 


3224 


> 24.5 


0.00 ± 0.20 


B 


VLT/FORS2 


4017 


> 25.9 


-0.08 ±0.08 


R 


VLT/FORS2 


5144 


> 23.6 


0.27 ±0.34 


z 


VLT/FORS2 


8135 


> 25.7 


0.02 ± 0.06 


i' 


Gemini-N/GMOS 


9350 


> 24.5 


0.02 ±0.18 


z' 


Gemini-N/GMOS 


10611 


22.80 ±0.16 


2.82 ± 0.44 


J 


Gemini-N/NIRI 


11785 


21.41 ±0.05 


10.21 ±0.50 


H 


Gemini-N/NIRI 


13280 


21.12 ±0.04 


13.26 ±0.51 


K 


Gemini-N/NIRI 


95658 


22.42 ±0.16 


4.0 ±0.6 


K 


Gemini-N/NIRI 


1.2 x 10 6 


> 27.07 




r' 


Gemini-N/GMOS 



Note. - Optical/NIR observations of GRB090429B. Magnitudes are 
quoted in the AB system, and corrected for the expected Galactic extinc- 
tion along the line of sight, Eb-v = 0.015. Quoted errors are la and limits 
are at the 3a level. 
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Table 3: Log of HST Observations of the GRB 090429B Field 



Date 


Start Time 


Inst/Filter 


Exp time 


Limit 


Flux density (/iJy) 


3 Jan 2010 


03:13 


ACS/F606W 


2100 


> 27.6 


0.005 ± 0.008 


10 Jan 2010 


21:54 


WFC3/F160W 


2412 






22 Feb 2010 


19:22 


WFC3/F160W 


2412 


> 27.5 


0.007 ± 0.005 


24 Feb 2010 


03:19 


WFC3/F105W 


2412 






28 Feb 2010 


13:56 


WFC3/F105W 


2412 


> 28.3 


-0.001 ± 0.005 



Note. — A log of the HST optical and NIR observations of the GRB 090429B field. Flux densities are given in 
the measured apertures and are not corrected for light outside the apertures. Errors are la and the limits are given 
in the AB- magnitude system at the 3a level (and do include aperture corrections). 



Table 4. Secondary Standards within the NIRI Field of View 



Star 


RA 


DEC 


^2MASS 


Jcai 


H-2MASS 


H ca / 


K-2MASS 


K ca z 


A 


14:02:35.05 


32:11:07.2 


14.754 


14.753 


14.411 


14.407 


14.212 


14.341 








±0.034 


±0.006 


±0.055 


±0.008 


±0.075 


±0.010 


B 


14:02:40.60 


32:09:28.9 


15.419 


15.451 


14.968 


15.023 


14.846 


14.944 








±0.055 


±0.007 


±0.076 


±0.009 


±0.127 


±0.015 


C 


14:02:38.11 


32:10:08.9 




19.523 




19.001 




18.585 










±0.035 




±0.028 




±0.032 



Note. - Vega magnitudes for our two secondary standard stars utilized in photometry 
of the afterglow of GRB 090429B. The 2MASS entries refer to the magnitudes contained 
within the 2MASS catalog, while those denoted cal refer to our improved values based on the 
UKIRT / WFC AM and ES02.2/GROND observations. Uncertainties are la. 
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Fig. 1. — X-ray (top) and optical/IR (bottom) lightcurve of GRB 090429B, the left-hand and 
bottom axis represent the observed time and flux/magnitude, while the top and right-hand axis 
show rest-frame time and luminosity, respectively. The solid points in the top panel show the 
observed XRT data, along with a solid line represent ing the model. The dashed line represents 
the best fit model for GRB 090423 (jTanvir et al.ll2009i ) overplotted as it would appear at z ~ 9.4 . 
The lower panel shows the optical lightcurve, along with a single power-law fit to the (red) i^-band 
points. (H and J are shown as green and blue, respectively. For clarity we have shown only the i- 
and z-band limits (cyan) in the optical). Additionally, the dashed line again shows the model of 
GRB 090423 at z ~ 9.4 . As can be seen, the luminosity and general behavior of GRB 090429B in 
both X-ray and optical is similar to that of GRB 090423. 
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Fig. 2. — Discovery images of the GRB 090429B afterglow. The images are all obtained from 
Gemini-North, and show the deep non-detection in the z band (which agrees with similar observa- 
tions in griz obtained at GROND, B, R, z obtained at the VLT, and an z-band image at Gemini), 
coupled with the relatively bright object seen in H and K. At z ~ 9.4 , Ly-a lies within the J 
band, and explains the marginal detection at that wavelength. 




Fig. 3. — IR to X-ray spectral energy distribution at To + 10 4 s can be explained by an intrinsic 
broken power-law spectrum. The green solid line extends the IR spectral slope derived from the 
fit to the optical/NIR data, albeit that the prior on /?o does essentially fix this value. The blue 
dot-dashed line extrapolates the unabsorbed X-ray spectrum to lower frequencies, showing that a 
single power-law fails to fit the broadband SED at this time. The red dashed line shows the SED for 
the best-fit extreme-redshift (z ~ 9.4 ) model, z' and i' upper limits are shown as black triangles. 



Fig. 4. — Wide-field image of the GRB 090429B field, obtained with Gemini-North 14 days after 
the burst. The location of the GRB is marked with a crosshair. Additionally, we mark the positions 
of the three comparison stars used to refine our photometry (note that star C is faint, and lies at 
the end of the marked arrow, just to the south of the galaxy), and the location of a large elliptical 
galaxy (Gl), which is the central galaxy of a modest cluster at z ~ 0.08, which may provide a 
modest lensing magnification. Note the silhouette of the guide probe obscures part of the field. 
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Fig. 5. — Input priors adopted for our photometric redshift fitting. [Left panel:] In the relativistic 
fireball model, the intrinsic spectral slope in the optical should lie between fix and Px — 0.5 (plus 
the associated measurement errors). To achieve this we use a lognormal distribution centered at 
0.5 (since there does appear to be a break between the optical and X-ray, see Figure [3]). This 
is a relatively weak prior and simply avoids extreme values of j3. Right panel: The second prior 
is on the intrinsic optical afterglow luminosity, and impacts solutions that would result in an 
unreasonably bright luminosity (it is not bounded at the faint end, and hence the low-redshift 
solutions are unaffected). It is therefore based on the empirically observed upper envelope of 
afterglow luminosities. The primary impact of this prior is to disfavor moderate {Ay > 3) scenarios 
at high redshift (z > 7), where the burst would have been more luminous than any other known 
afterglow. 



-27- 




Fig. 6. — Spectral energy distribution of the GRB 090429B afterglow formed by extrapolating our 
observed photometry to 3 hr post-burst assuming the magnitude remains constant, i.e. a = (for 
varying a fits see Figure [8|). The vertical error bars represent la uncertainty, and the horizontal 
shaded bars illustrate the widths of the broadband filters. The best fit model (x 2 /dof = 1.76/3) 
to the data points is shown as the solid red line, the parameters being redshift z = 9.36, rest- 
frame extinction Ay = 0.10 and intrinsic power-law slope /3o = 0.51. The inset simply replots the 
short wavelength part of the figure (indicated by a dotted box) on a logarithmic flux density scale, 
to more clearly show the constraints from the optical measurements. An alternative low-redshift 
(z ~ 0), high extinction (Ay = 10.6) model is shown as a dashed blue line, but in fact is formally 
ruled out at high significance (x 2 /dof = 26.2/4). 
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a= 0.00 
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Fig. 7. — Confidence contours on a parameter space of redshift and host galaxy extinction for the 
GRB 090429B afterglow, for our favored set of prior assumptions (green contours are 90%, 99% and 
99.9% confidence). The gray scale shows likelihood down to much lower levels, formally ~ 10 -7 . 
All fits at z < 7.7 are ruled out at > 99% confidence, and while fits can be found at z ~ they 
are markedly worse than the high-z solutions. The best z < 5 solution (formally at z = 0) is 
marked with the blue cross and requires Ay ~ 10, and is also disfavored by the lack of any host 
galaxy to deep limits, and the inconsistency of the required Ay with the hydrogen column density 
measured from the X-ray afterglow. To illustrate this the best-fit A^h from the X-ray spectrum is 
converted into Ay and plotted onto the contour plot as the purple lines (dashed lines show the 90% 
error range, and the dotted li nes show the limits of the systematic error due spanning the range of 
gas-to-dust ratios reported bv lSchadv et al.ll2010l ). As can be seen the Ay inferred from the X-ray, 
and that required from the photometric redshift fit are inconsistent at low redshift, but broadly 
consistent with the high-z fit. 
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Fig. 8. — Results of SED fits with a range of possible values for the temporal power-law index, and 
different reddening laws. The plots encompass the canonical reddening laws for the SMC, LMC and 
Milky Way (which are chara cterized by the increasing influence of the 2175 A bump) as well as the 
law of lMaiolino et al.l (|2004l ) which is approximately flat ("gray") from ~ 1800 to 3000 A. As can 
be seen, the assumed temporal index has only a minimal impact on our results, and our assumption 
of a = therefore does not affect our analysis. The majority of GRB aftergl ows are best fit wit h 
SMC-like absorption, and we therefore adopt this as our choice model (e.g. ISchadv et al.l 120071 ). 
Other laws can produce broader allowed redshift ranges, in particular extending as low as z ~ 6.3 
at 99% confidence in the lMaiolino et al.l (|2004l ) case, but all rule out low-z (< 6) scenarios. 
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Fig. 9. — Posterior likelihood plotted on both a linear (upper) and log (lower) scale, for the models 
assuming an SMC dust law, where we have marginalized over both a (assumed to have a flat prior 
between —1 and +1) and Ay (assumed to have a flat prior between and 12). The dark and light 
shaded bands show the extent of the 90% and 99% enclosed likelihood regions respectively. 
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Fig. 10.— Our late time HST observations of the GRB 090429B field in the optical and NIR. No 
host galaxy is detected in any filter, supporting a high-redshift origin, since a host with z < 1 
would be very unlikely to be fainter than these limits, even if dusty. At F160W the host remains 
undetected, but the observations reach limits which would uncover ~ 50% of the z > 8 candidates 
in the Hubble Ultra-Deep Field (UDF). Hence, the non-detection of any host is fully consistent 
with our high-z model, but inconsistent with any lower redshift, high extinction scenario. 
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Fig. 11. — Solid lines show the 3a absolute magnitude limits for GRB 090429B in each of our filters 
F606W (blue), F105W (gree n) and F160W (red) . The inferred absolute magnitudes (AB) of a 
sample of GRB h ost galaxies (jFruchter et al.ll2006l ). as a function of redshift (plot modified from 
Perlev et al.ll2009i ). The known GRB hosts with HST observations are p lotted as red points , and 
are supplemented at high redshift by the observations of GRB 050904 bv lBerger et al.l ([20061 ). As 
can be seen, all of these lines lie significantly below the majority of GRB hosts and offer support 
for a high -redshift origin f or GR B 090429B. The blue points at z ~ 8 represents the Lyman break 
sample of iBouwens et al.l (120101 ) . As can be seen, the limiting magnitude for GRB 090429B lies 
roughly at the median of this distribution, and so the non-detection in our observations would not 
be unexpected at z ~ 9.4 



